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ABSTRACT 

In this paper we study the properties of pseudobulges (bulges that appear similar to disk galaxies) and clas- 
sical bulges (bulges which appear similar to E-type galaxies) in bulge-disk decompositions. We show that the 
distribution of bulge Sersic indices, nb, is bimodal, and this bimodality correlates with the morphology of the 
bulge. Pseudobulges have «/> < 2 and classical bulges have «6 > 2 with little-to-no overlap. Also, pseudobulges 
do not follow the correlations of Sersic index with structural parameters or the photometric projections of the 
fundamental plane in the same way as classical bulges and elliptical galaxies do. We find that pseudobulges 
are systematically flatter than classical bulges and thus more disk-like in both their morphology and shape. 
We do not find significant differences between different bulge morphologies that we are collectively calling 
pseudobulges (nuclear spirals, nuclear rings, nuclear bars, and nuclear patchiness) appear to behave similarly 
in all parameter correlations. In Sersic index, flattening, and bulge-to-total ratio, the distinction appears to be 
between classical bulges and pseudobulges, not between different pseudobulge morphologies. The Sersic in- 
dex of the pseudobulges does not correlate with B/T, in contrast to classical bulges. Also, the half-light radius 
of the pseudobulge correlates with the scale length of the disk; this is not the case for classical bulges. The 
correlation of Sersic index and scale lengths with bulge morphology suggests that secular evolution is creat- 
ing pseudobulges with low-Sersic index, and that other processes (e.g. major mergers) are responsible for the 
higher Sersic index in classical bulges and elliptical galaxies. 

Subject headings: galaxies: bulges — galaxies: formation — galaxies: structure — galaxies: fundamental 
parameters 



L INTRODUCTION 

Historically, all bulges were assumed to be little elliptical 
galaxies residing in the centers of galactic disks. Thus it was 
also assumed that all bulges were dynamically hot stellar sys- 
tems. Recent work puts this assumption in question. Many 
bulges have disk-like features that do not resemble E-type 
galaxies. 

Classical bulges are dynamically hot, and relatively fea- 
tureless; they appear similar to the end products of galaxy 
major mergers. They are easily recognized as having mor- 
phologies very similar to E-type galaxies. 

In contrast, many bulges have features that are not found 
in elliptical galaxies, but in galactic disks. These features 
include the following: kinematics dominated by rotation 
(Kormendv 1993); flatte ning similar to th at of their outer disk 
ffathi & Peletier ,20031 lKormendvllI993i) ; mi d-IR colors o f 
the bulge ar e similar to those of the outer disk (lFishedl2006h : 
nuclear bar (Erwin & Sparke"2002); nuclear ring and/or nu- 
clear spiral (Carollo et al. 1997); nea r-exp onential surface 
brightness profiles ( Andredakis & Sanders! [l994l) . Bulges 
with such properties are called pseudobulges. All of these 
phenomena are manifestations of stellar systems that that are 
dynamically cold. However, the extent to which these fea- 
tures exist simultaneously in all pseudobulges is not yet well- 
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un derstood. For a review of the p roperties of pseudobulges 
see lKormendv & Kennicutll (|2004|) (KK04 hereafter). 

In this paper, we refer to the super-set of the two systems 
simply as bulges. We note that in this paper the term pseu- 
dobulge refers to bulges with morphology reminiscent of disk 
galaxies; there is no a priori assumption about their formation 
mechanism. 

Investigating the link between bulge morphology and bulge 
structural parameters is the main interest of this article. More 
precisely, those bulges with morphologies that are reminis- 
cent of disks (such as nuclear rings, nuclear bars and nu- 
clear spirals) are expected to have lower Sersic index than 
those bulges with smooth light distributions resembling E- 
type galaxies. We will investigate whether this expectation 
holds and whether the distribution of bulge Sersic indices is 
dichotomous. 

Determining whether a critical Sersic index that discrimi- 
nates between classical bulges and pseudobulges exists, and 
what its value is, impacts our understanding of bulge for- 
mation in at least two ways. First, the distribution of Ser- 
sic indices would constrain formation theories of classical 
bulges and pseudobulges. Also, the existence of a critical Ser- 
sic index would robustly establish a method for distinguish- 
ing pseudobulges from classical bulges without using high- 
resolution imaging or kinematic data, neither of which is cur- 
rently available in large surveys. 

In the last ten years, the subject of surface brightness pro- 
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files of bulges and elliptical galaxies has experienced a shift of 
paradigm. Traditionally, surface brightness profiles of ellipti- 
cal galaxies and bulges were all thought to be of a single shape 
that is well-chara cterized by the de Vaucauleurs r'/"^ profile. 
ICaon et aTl (1 19941) show that surface brightness profiles of el - 
liptical galaxies are better fit by Sersic profiles (ISersidI 19681) . 
which generalize the ex ponent in the de Vaucauleurs p rofile 
to a free parameter Also. lAndredakis & Sandersl(ll994!) show 
that many bulge-disk galaxies are better described by a dou- 
ble exp onential than an inner r^l ^ profile with an outer expo- 
nential. 'Andredakis et al.l (11995 ^ generalize this to show that 
bulge-disk galaxies are better fit by inner Sersic profiles with 
outer exponential disks than double exponentials profiles. The 
Sersic function plus outer disk model for bulge-disk galaxies 
reads 



/(r) = /()exp 



(r/ro)'/""l+/rfexp[-(r//z)] , (1) 



where /o and ro represent the central surface brightness and 
scale length of the bulge, /,/ and h represent the central surface 
brightness and scale length of the outer disk, and rih represents 
the bulge Sersic index. 

It is not surprising that Sersic profiles fit bulge surface 
brightness profiles better than r'/^'-profiles, since the Ser- 
sic function has more flexibility due to the extra parame- 
ter However, the new parameter, Ub, correlates with many 
properties of the stellar systems to which it is fit, includ- 
ing but not limited to the following: velocity dispersion a 
jKhosroshahi et al. 2000) , absolute magnitude (Graham et al. 
[19961) . and effective radius (ICaon et al.l 119941) . Many au- 
thors have shown that these c orrelat i ons extend to bulges of 
bulge-disk galaxies (Graham llOOlt iMacArfliur et al.1 1200I 
de Jong et al. '2004'; Tho mas & DaviesI 12006"). Additionallv. 
Andredakis et al. (1995) show that the Sersic index of bulges 
correlates with Hubble type, decreasing from rib ^ 3.1 for SO 
galaxies to nb ^ 1.6 in Sbc-Sd galaxies. Therefore it is rea- 
sonable to assume that it has physical significance. The Sersic 
index is often referred to as the shape parameter, as it is gener- 
ally taken as a surrogate for properties such as concentration 
of the surface brightn ess profile. For a review o f the properties 
of Sersic profiles see lGraham & Driver! (120051) . 

The tentative assumption is that the Sersic index of a bulge 
reflects the classical bulge - pseudobulge dichotomy. Lower 
Sersic index might indicate that a bulge is more likely to be a 
pseudobulge. We do not understand the mechanism that is re- 
sponsible for determining the Sersic indices in pseudobulges 
(or classical bulges). Yet, it seems plausible that the light dis- 
tribution be similar to that of a disk, since so many other of 
the propert ies of pseu dobulges are similar to those of galac- 
tic disks. Courteau et al.l (1199 6) use bulge-disk decomposi- 
tions of 243 galaxies to show that the 85% of bulges in Sb-Sc 
galaxies are better fit by the double exponential than cuspier 
r^^^ models. Thus, the common conclusion is that pseudob- 
ulges are mar ked by near exponential Sersic index (KK04). 
IScarlata et al.l (12004) shows with STIS acquisition images that 
bulges with surface brightness profiles more resembling expo- 
nential profiles are more likely to have disk-like morphology 
(e.g. spiral arms), yet there is significant scatter to this claim, 
in their sample. They go on to show that the distribution of 
central slopes of surface brightness profiles of the bulges in 
their sample is bimodal when plotted against absolute magni- 
tude. 

It is yet unknown where pseudobulges should lie in other 
parameter correlations, such as fundamental plane projec- 



tions. We do not expect that pseudobulges occupy a signif- 
icantly different location than classical bulges in fundamen- 
tal plane parameter space, since the fundamental plane is not 
known to be bimodal. Many studies of th e locations of bulges 



in structural p a rameter space ex i st (e.g . iBender et al.l 
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Graham' "200?, 'MacArth ur et"an llOOl Ide Jong et al.l l2004t 
Thomas & Davies 2006), and no significant bimodal behavior 
is noticed. Carollo ( 1999) remarks, though, that pseudobulges 
deviate more from the ^e-re relation (Kormendv 1977 ). We 
will investigate this further in this paper. Kormendv (1993h 
shows that the majority of bulges are rounder than the outer 
parts of the disks they reside in, yet a significant minority are 
as flat as their associated outer disk. A few bulges are even 
flatter than their outer disk. This behavior correlates with 
Hubble type; bulges in later type galaxies hav e flattening more 
similar to that of the associated outer disk. i Fathi & Peletied 
([2003) carried out bulge-disk decompositions of 70 galaxies 
on higher-resolution data finding a similar result. 

As discussed in KK04, pseudobulges are characterized 
principally by having less random motion per unit stellar 
light. They are rotation-domin ated systems ( Kormendvll993l: 
iKormendv & Kennicuttll2004h . Thus it makes sense that they 
be flatter The relative flatness of a bulge to its associated 
outer disk has also been suggested as a pseudobulge indicator 
(KK04). We will test this hypothesis in this paper by compar- 
ing the flatness of bulges with disk like morphologies to that 
of bulges with morphologies like those of elliptical galaxies. 

The properties of pseudobulges are in stark contrast to the 
expected end result of the hierarchical merging process; one 
does not expect violent relaxation to produce spiral structure 
and dynamics that is dominated by ordered motion. Fur- 
ther scenarions for the formation of bulges have been sug- 
gested. Clump instabilities in disks at high redshift can form 
bulge-like structures in simulations (Noguchi 1999). It is also 
plausible that gas rich accretion could form dynamically cold 
bulges. Internal evolution of disks can drive gas and stars to 
the center of a disk galaxy as well. The population of bulges 
as a whole and any one particular bulge may be the result of 
more than one of these processes. 

However, the connections between bulge and disk stel- 
lar po pulations (iPeletier & BalceU sI Il 99(5t iMacArthur et all 
120041) . inter-steUar medium dRe gan et alT |2001l; lFishedl2006l) 
and scale lengths (^ Courteauet al.i il996h may suggest that 
pseudobulges form through processes intim ately linked to 
their host disks. Furthermore. iDrorv & Fisher (i2007h find that 
classical bulges occur in red-sequenc e galaxies and pseudob- 
ulges occur in blue cloud galaxies. iKormendy & Kennicut j 
('2004|) reviews the case that pseudobulges are not the result 
of major mergers, but rather that internal disk-evolution may 
be responsible for them. However, that bulges can form out 
of the int ernal evolut ion of a disk is not a new idea; see, for 
example,|Hohl(|T975|). 

Simulations suggest that if disk galaxies do not experi- 
ence major mergers, they may evolve by redistribution of en- 
ergy and angular momentum driven by non-axisymmetries 
such as bars, ovals, and spiral structure (ISimkin et al.] 1 19801; 



iPfenniger & Friedlilll99l1:lDebattista et alJl2004 . resulting in 
star formation and bulge-like stellar densities, thus forming 
pseudobulges. Indeed, a correlation between central star for- 
mation rate and the p resence of bars and ovals has been de- 
tected (e^g^piethetal. 2005; Jo gee et a l. 2005; Fisher 2006J. 
Also, iPeepl es & Martini ( 200^ find that galaxies with nu- 
clear rings and/or nuclear spirals are more strongly barred. 
Thus we also investigate the possible connection between 
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driving mechanisms and structural properties. 

It is, thus, possible that the absence of a classical bulge in 
a galaxy indicates that the galaxy has not experienced a ma- 
jor merger since the formation of the disk. In this context, 
pseudobulges may be thought of as more similar to pure disk 
galaxies that have a surface brightness profile which breaks 
from the outer exponential profile to a more steep inner sur- 
face brightness profile. Though the frequency of pseudob- 
ulges has not yet been robustly calculated, if they are common 
then this implies that many disk galaxies did not suffer major 
mergers since their formation. 

This paper is organized as follows. In § 2 we present the 
observational data we use, and we discuss the surface bright- 
ness fitting procedure. In § 3 we present results on the loca- 
tion of pseudobulges and classical bulges in various structural 
parameter correlations. In § 4 we discuss the flatness of pseu- 
dobulges and classical bulges. In § 5, we discuss behavior 
of different bulge morphologies (nuclear bars, nuclear spirals, 
and nuclear rings) in various parameter correlations. In § 6 we 
summarize and discuss these results. Finally, the appendix 
includes an image of each galaxy, all decompositions and a 
discussion of the robustness of our decomposition and fitting 
procedure. 

2. METHODS AND OBSERVATIONS 
2.1. The Sample 

The aim of this work is to establish whether or not pseu- 
dobulges - recognized by the presence of disk-like morpho- 
logical features as motivated and discussed in KK04 - can 
be distinguished from classical bulges simply by structural 
features in their surface brightness profiles, most prominently 
their profile shape. Thus we ask whether bulges that contain 
disk-like morphologies (pseudobulges) have lower Sersic in- 
dex and higher flattening ratios than bulges with elliptical- 
like morphologies (classical bulges). Answering this ques- 
tion requires high resolution imaging (preferably in the opti- 
cal bands) to detect the nuclear spirals, bars, and rings; and 
we need surface brightness profiles with large dynamic range 
in radius to accurately determine the parameters in Eq. 1 for a 
bulge-disk decomposition. 

We choose galaxies observable from the northern hemi- 
sphere that have data in the Hubble Space Telescope (HST) 
archive. We limit our selection to galaxies closer than ~ 
40 Mpc to resolve features in the bulge. 

The link between non-axisymmetries (barred and oval dis- 
tortions) and secular evolution motivates us to create a sam- 
ple containing roughly equal numbers of galaxies with a driv- 
ing agent (galaxies with a bar and/or an oval) and galaxies 
without a driving age n t. Det ection of oval distortions are dis- 
cussed in iKormendyl (Il982h . They are identified by nested 
shelves in the surface brightness profile usually having dif- 
ferent position angles. We identify bars by consulting the 
Carn egie Atlas of Galaxies (Sandage & Bedke 1994) and the 
RC3 (Ide Vaucouleurs et al.lll99 1). If a galaxy has both a bar 
and an oval, we call that galaxy barred. Additionally, we look 
for bars and ovals in all galaxi es using K-band images from 
2MASS (S krutskie et al.ll2006l) . Note that we do not distin- 
guish grand design spirals as a possible secular driver, though 
they may be able to generate a similar but less extreme effect 
as bars do (KK04). We use 39 undriven (no bar and no oval) 
galaxies and 38 driven galaxies (30 barred and 8 ovaled), a 
total of 77 galaxies. 

In Fig. [U we show the distribution of global properties of the 



galaxies in our sample; these are also listed in Table 1 . The 
distribution of the distances of the galaxies in our sample is 
heavily peaked at 16 Mpc due to the Virgo cluster and has a 
standard deviation of 6 Mpc. We derive total magnitudes by 
2D integration of our surface brightness profiles. The distri- 
bution of absolute V magnitudes ranges mostly from -19 to 
-22 with a median value of -20.5. 

KK04 compiles data from several different studies to gen- 
erate preliminary statistics on how the frequency of pseudob- 
ulges varies along the Hubble Sequence. They suggest that 
both pseudobulges and classical bulges exist at intermediate 
Hubble Types (SO to Sbc). There appears to be a transition 
from classical bulges being more frequent at early types (SO- 
Sb) to pseudobulges being more frequent at later types (Sbc - 
Sd). They further suggest that classical bulges will be almost 
non-existent at Hubble types Sc and later. In the right panel 
of Fig. [H we show the distribution of Hubble Types (taken 
from Sandage & Bedkel ll994<) of the galaxies in our sample. 
To test for differences between pseudobulges and classical 
bulges we choose to sample the range of Hubble type from 
SO to Sc. If we combine this with our choice of evenly sam- 
pling driven and undriven galaxies, we expect that our sam- 
ple should overemphasize pseudobulges. The distribution of 
Hubble types is as foUows: 13 SO, 2 SOa, 18 Sa, 6 Sab, 20 
Sb, II Sbc, 9 Sc. We will also use 24 E type galaxies from 
iKormendy et al.l (12006) as early-type sample in some param- 
eter correlations. 

2.2. Identification of pseudobulges 

We define "bulges" photometrically, as excess light over the 
inward extrapolation of the surface brightness profile of the 
outer disk. The region of the galaxy where this excess light 
dominates the profile is the bulge region. We classify galax- 
ies as having a pseudobulge by their morphology within this 
bulge region; if the bulge is or contains any of the follow- 
ing features: a nuclear bar, a nuclear spiral, and/or a nuclear 
ring, then the bulge is called a pseudobulge. Conversely, if 
the bulge better resembles an elliptical galaxy (relatively fea- 
tureless isophotes), then the bulge is called a classical bulge. 
This method is discussed in KK04. The existence/absence 
of visibly identifiable disk-like structure in a bulge correlates 
with properties of the bulge and the whole galaxy. The same 
method is shown to be successful in identifying bulges with 
higher specific star formation rates (F isher 2006) and globally 
bluer galaxies (iDrory & Fisheijl20()7l) . 

Fig. |2] shows high-resolution HST images of the bulge re- 
gion of all galaxies in our sample. All images are taken in 
close-to V band filters. Note first that classical bulges (for ex- 
ample NGC 1398, NGC 2775, NGC 2880, and NGC 3115, all 
in the rightmost column of the first page of the figure) have a 
smooth stellar light profile. There is no reason evident in the 
images to think that any of these galaxies harbor a pseudob- 
ulge. The bulges fit the description of E-type galaxies. On the 
other hand, NGC 4030 shows a face-on nuclear spiral (bot- 
tom row on the second page). The spiral dominates the radial 
profile for more than a kiloparsec. NGC 4736 (fourth page, 
second row) also has a nuclear spiral pattern but also has a 
nuclear bar; note how the spiral arms seem pinched vertically 
in the image. Moellenhoff et al. ( 1995) study this nuclear bar 
in more detail using dynamical modeling. NGC 4371 (third 
page, middle row) is an SO galaxy with at least one nuclear 
bar. NGC 3351 (second page, third row) has a prominent nu- 
clear ring that heavily distorts the surface brightness profile 
(see Fig. |4li; this nuclear ring is quite well known (.Sandag3 
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Fig. 1 . — The distribution of absolute magnitude (left), distance modulus (middle), and Hubble Type (right) for all 77 galaxies in our sample. 



[196 IV The bulge of NGC 2903 (first page, bottom row) shows 
a chaotic nuclear region; it appears nearly spiral but is not reg- 
ular enough over a significant radial range to call it a nuclear 
spiral. 

In this paper we classify bulges with near-V-band images 
(F547M, F555W, & F606W). Thus our method is subject to 
the effects of dust obscuration. However, the structures used 
to identify pseudobu lges are usua lly experiencing enhanced 
star formation rates (lFisheiil200 6l). Pseudobulges are, there- 
fore, easier to detect in the optical band passes where the 
mass-to-light ratios are more affected by young stellar pop- 
ulations, rather than in the near infrared where the effects of 
dust are less pronounced. It is important to note that classi- 
cal bulge s may have dust in their center, as do many elliptical 
galaxies jLauer et al.l l2005h . In fact, many classical bulges 
shown in Fig. |2]have some dust lanes in their bulges. The 
presence of dust alone is not enough to classify a galaxy as 
containing a pseudobulge; instead, it must be of a disk-like 
nature. 

2.3. Data Sources and Surface Photometry 

As stated in the introduction, the Sersic function provides 
better fits and more information over two-parameter fitting 
functions describing surface brightness profiles of bulges and 
elliptical galaxies. However, this infor mation come s at the 
price of more detailed observations. ISaglia et"an ( Il997h 
shows that replacing an r'/'^-model with a Sersic model 
makes little difference in residuals, unless one has data with 
high dynamic range in radius. Further, the coupling be- 
tween parameters in the Sersic function can be quite high 
jGraham & Collesslll997h . Thus again, it is necessary to fit 
the decomposition (Eq. 1) to large radial range in order to 
minimize these degeneracies. For each galaxy we therefore 
combine multiple data sources together: high-resolution HST 
imaging in the center, and ground based wide-field images 
covering the outer disk. Comments on data sources follow. 
Table 1 lists the sources of data used for each galaxy. 

All profiles contain HST data sources. PC2 data has a 
small field-of-view ('^ 18 x 18 arcsec^); thus, it is critical to 
supplement PC2 data with wide field data. ACS/WFC has 
proven to be an excellent instrument for obtaining large ra- 
dial fitting range. It provides a reasonable sized field-of-view 
(~ 100 X 100 arcsec^) at high spatial resolution (0.049 arcsec 



pixel"'). 

For as many galaxies as possible we obtain wide field im- 
ages from the Prime-Focus-Camera on the McDonald 0.8 m 
telescope. This instrument provides a large unvignetted field 
of view (45 x 45 arcmin^), and a single CCD detector There- 
fore, we can more accurately carry out sky-subtraction. These 
data generally are the deepest. We also include images 
from the Sloan Digital Sky Survey ( Ade lman-McCarthv et alj 
12006) . 2Micron All Sky Survey (2MASS) and the Isaac New- 
ton Group (ING) Archive. Individual data sources are noted 
in Table 1 . 

All raw data (McDonald 0.8 m & ING data) ai-e bias- 
subtracted, flat-fielded, and illumination corrected. We sub- 
tract the sky background by fitting a plane to a smoothed ver- 
sion of the image where the galaxy and bright stars have been 
removed. 

We calculate Johnson V-band magn itude zero points using 
the transfor mations in Holtzma n et alj (|l995) for the WFPC2 
images and ISirianni et al.. (2005F for the ACS images. SDSS 
g and r profiles are converted to a singl e V-band profile fo r 
each galaxy using the transformations in lSmith et al.l (120021) . 
Other profiles are simply shifted to match in surface bright- 
ness. These transformations are derived on galactic disk stars, 
not galaxies. Further, the calculations require color informa- 
tion. We use colors from Hyper-LEDA, which refer to colors 
of the entire galaxies, and the galaxies in our sample most cer- 
tainly have non-zero color gradients. Therefore the absolute 
values of surface brightness in this paper are not expected to 
be consistent to more than 0.3 mag. However, this does not 
affect our conclusions which are based the structure in the 
profiles and not on absolute magnitude. We check that our to- 
tal magnitudes are consistent with those published in the RC3 
and Hyper-LEDA. 

We use the isopho te fitting routine of 
iBender & MoeUenhoffI (Il987h . We identify and mask 
interfering foreground objects in each image. Then we fit 
ellipses to isophotes by least squares. Here, isophotes are 
sampled by 256 points equally spaced in an angle 6 relating 
to polar angle by tan 9 = (a /b) tan (p, where is the polar 
angle and b/a is the axial ratio. The software determines 
six parameters for each ellipse: relative surface brightness, 
center position, major and minor axis lengths, and position 
angle along the major axis. 
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Fig. 2. — High-resolution HST images of all bulges in our sample. The images are scaled logarithmically, and the range of shown intensities is chosen to 
emphasize those features which motivate the bulge classification. For all galaxies we give the galaxy name, imaging instrument used, filter, and the classification. 
We label bulges classified as classical by 'c' and pseudobulges by 'p'. The images show the inner 40 arcsec (ACS), 36 arcsec (PC2), and 80 arcsec (WF2) of the 
galaxies. 
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Fig. 3. — Comparison of structural parameters (ro&nj,) derived from bulge- 
disk decompositions in the //-band to those in the V-band. The identity rela- 
tion is marked by the dashed line. 



2.4. Bulge-Disk Decomposition 

We carry out a bulge-disk decomposition on each galaxy 
in our sample by fitting Eq. 1 to the major axis surface 
brightness profiles. Ou r decomposition code is also used by 
iKormendv et alj ( l2006l) to fit SO and E-type galaxies with the 
disk component "turned off". The average root-mean-square 
deviation is ~ 0.09 ±0.03 mag arcsec"^. The largest deviation 
of any of the data from its fitted profile is 0.18 mag arcsec"^. 

The parameters determined in bulge-dis k decomposition 
with a Sersic model bulge are coupled. iMac Arthur et alj 
( l2003h shows that if the bulge is sufficiently large and the 
resolution is sufficiently small, initial parameter estimates do 
not affect the final fit too much. For the most part, this is 
true. Yet, our experience is that initial parameter estimates 
may still affect the resulting fit. For each galaxy, we begin 
with a large parameter range that is symmetric about the ini- 
tial guess. Then we re-fit the galaxy iteratively adjusting the 
range of allowed parameters to be narrower. This results in 
slight changes of best-fit values with lower than without 
this iteration. The details of each profile are used to decide the 
width of available parameter space each parameter is given. 
Typically the available range for rtt, is Ant, ~ 2-3. 



The decomposition is carried out on a major axis profile us- 
ing the mean isophote brightness. It does not take ellipticity 
into account during the fitting. However, these galaxi es are 
known to have varyin g ellipticity profiles CKormendvil 19931 : 
lFathi&Peletie3l2003h . Thus, we take the mean ellipticity 
for each component and adjust the luminosity accordingly: 
L = (1 -e)Lfif The radius of the component is defined as the 
radius range within which that component dominates the light 
of the profile. We also adjust all magnitude dependent quan- 
tities (luminosit y and surface br i ghtn ess) for Galactic redden- 
ing ac cording tolSchlegel et al.l (1199 8). We use the distances 



of Tull v~& Fishend 19881) . We do not make any corrections for 
extinction within the galaxies being studied. 

We carry out the bulge-disk decompositions in the V band. 
This ensures that our results will remain applicable to large 
surveys commonly done in the optical bands (e.g. SDSS). 
However, the radial variation in age and metallicity that ex- 
ists within gala xies may bias the parameters of bulge-disk 
decomposition. ' Carollo et alj (2002) shows that the centers 
of intermediate type galaxies contain significant structure in 
V -I color maps, and this variation occurs on scales smaller 
than the bulge. Thus, there is doubt as to whether the parame- 
ters derived on profiles in the middle optical band passes truly 
reflect the properties of the stellar mass distribution. 

It is beyond the scope of this paper to investigate the cor- 
relations of colors with Sersic parameters. Yet, as a check 
of the stability of structural quantities compared to color gra- 
dients, we carry out bulge-disk decompositions on all galax- 
ies in our sample that have archival NICMOS images in the 
filter F16QW {H band). We supplement those data with 
ground-based wide-field data. For this purpose we use mainly 
2MASS data and any H-band data available in NED. We then 
compare the fit parameters in V to those in H. The results are 
shown in Fig. [3] In those parameters which do not depend 
on magnitude (ro and «;,) there is little difference. The av- 
erage difference is A,, ^n;, ^ -0.03, and A,,^ log(ro) ~ 0.07. 
This is similar to the results of other papers that have done 
decompositions at multiple wavelengths (e.g. lGrahamll200il : 
iMacArthur et al. 2003). 

In Fig.|4]we show four examples of our photometry and the 
bulge-disk decompositions. For each galaxy we show an ellip- 
ticity profile, a residual profile, and a surface brightness pro- 
file along with the decomposition determined by fitting Eq. 1 
to the galaxy's surface brightness profile. These examples are 
selected to show the range in typical profiles in this paper. The 
entire sample of fits is shown in the appendix. 

The top two panels of Fig. |4] shows two galaxies 
(NGC 3031, left panel: NGC 3169, right panel) that are well 
described by an inner Sersic function and an outer exponen- 
tial disk. In both profiles the fit covers a dynamic range of 
~12 mag arcsec-2. NGC 3031 is well fit (RMS=0.06 mag 
arcsec"^) by Eq. 1 over Alog(r/arcsec) = 3.3 and NGC 3169 
for Alog(r/arcsec) = 2.51. Deviations are small, and there 
is little-to-no substructure evident in these profiles. Note that 
each bulge has a nuclear excess of light over the inward ex- 
trapolation of the fit, and that this nucleus is excluded from 
the fit (vertical hatches indicate beginning and end of fit range 
and excluded data points are marked by open symbols). 

Despite its successes, the Sersic bulge plus outer exponen- 
tial disk model of bulge-disk galaxies does not account for 
many features of galaxy surface brightness profiles. Disks 
of intermediate type galaxies commonly have features such 
as bars, rings, and lenses (see 'Kormendy"1982' for a descrip- 
tion of these). Further. Carollo et al. (.2002) shows that many 
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Fig. 4. — We show four example decompositions, spanning the range of the types of profiles that we find in the galaxies investigated in this paper. All data are 
plotted against r'/'' as this projection emphasizes deviations from the fit better than log(r), and shows the central regions better than plotting against Hnear radius 
would. For each galaxy we show a profile (black dots) of (from bottom to top panel) V band surface brightness; residual profile of the data minus -2.5 log /() ), 
where l{r) is the fit of the data to Eq. 1 ; and ellipticity profile, where c = l—h/a, and h/a is the axial ratio of the galaxy. In the bottom panel we also show the 
decomposition: the inner Sersic profile and the outer exponential profile as blaclt lines. The vertical hatches on the profile indicate the boundaries of the regions 
included in the fit, the bottom caption also notes the total range of the fit. Excluded data regions are marked by open symbols, fitted regions by closed symbols. 
The caption gives the galaxy name, the Hubble type, the total absolute magnitude, the Sersic index of the bulge and the root-mean-squared deviation of the fit to 
the fit region of the data. 
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bulges of early and intermediate type galaxies contain nu- 
clei. The bottom panels of Fig. |4] shows two such examples 
of galaxy profiles with significant deviations from the fitted 
decompositions. NGC 4448 (bottom left) is an example of 
a barred galaxy in which the bar is an especially prominent 
perturbation to the outer exponential surface brightness pro- 
file. NGC 3351 (bottom right) is an example of a compli- 
cated surface brightness profile, with multiple substructures 
that are not well described by the smooth nature of the bulge- 
disk model used here. This galaxy contains a nuclear ring 
near ~ 4 arcsec, and a bar from about 20 to 80 arcsec. These 
galaxies are not well described by Eq. 1, yet we do our best to 
decompose as many galaxies as possible. 

Bars, rings, lenses, and similar features do not conform to 
the smooth nature of Eq. 1, hence we carefully exclude re- 
gions of the profile perturbed by such structures from the fit. 
This is a risky procedure, as it requires selectively removing 
data from a galaxy's profile, and undoubtedly has an effect 
on the resulting parameters. We are often helped to identify 
bars by the structure of the ellipticity profile. Notice the the 
peak in ellipticity near 50 arcsec in NGC 3351. We try to 
err on the side of removing the fewest points possible. For 
those galaxies in which a bar is present, it is our assumption 
that removing the bar from the fit provides the best estima- 
tion of the properties of the underlying bulge and disk. We 
test our method by removing a typical number of points from 
a few galaxies with smooth profiles (NGC 2841, NGC 1425, 
NGC 4203, & NGC 1325). The result is that Sersic index 
can vary as much as An^ ~ 0.5, and characteristic radius 
by Alog(ro) ^ 0.5. The variance of these two parameters is 
tightly coupled and this is reflected by the uncertainty. Cen- 
tral surface brightness was mostly unaffected. If a region is 
not included in a fit we show that in the figure by using open 
symbols. A detailed discussion of the effects of bars and other 
features on the surface brightness profiles is given in the ap- 
pendix. 

3. STRUCTURAL PROPERTIES OF PSEUDOBULGES 
AND CLASSICAL BULGES 

3.1. Bulge Prominence Of Classical & Pseudobulges Along 
The Hubble Sequence 

The primary distinction between most previous studies and 
this one is that we do not treat bulges as an homogeneous set 
of objects. Here, we report on how morphologically identi- 
fied pseudobulges and classical bulges, taken to be distinct 
entities, behave in the parameters obtained from bulge-disk 
decomposition. Further we wish to know if the Sersic index 
is able to distinguish these two separate bulge types as has 
been suggested by many authors and not yet systematically 
tested. 

Fig |5] shows how the bulge-to-total ratio (luminosity of the 
bulge, divided by total luminosity of the fit; hereafter B/T), 
and the ratio of the bulge half-light radius (r^) to the scale- 
length of the outer disk (h) both correlate with Hubble type. 
We calculate the effective radius as 



: (fe«)'Vo, 



(2) 



w here b„ is a propo r tional ity constant whose expansion given 
in lMacArthur et all (12003). In all figures red circles represent 
classical bulges and blue crosses represent pseudobulges. In 
fig 12] the connected black squares show the average for each 
Hubble type. There are only two SO/a (T=0) galaxies; thus 
the dip in the average of B/T could merely be small number 
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Fig. 5. — The bulge-to-total ratio, B/T, top, and the ratio of the effec- 
tive radius to the disk scale length, re/h, bottom, of classical bulges (red 
circles) and pseudobulges (blue x's) plotted against Hubble types. Tlie black 
squares, connected by the black line, represents the average for all bulges 
(both pseudo- and classical) in each type. 

Statistics. Also, a few of the Sc (T=5) galaxies are not shown 
due to their very small B/T. 

As shown in the top panel of Fig.|5] bulge-to-total ratio is a 
decreasing function of Hubble type. This is not surprising as 
B/T is part of the original classification criteria. This behav- 
ior has been fo und by many authors, for a detailed study see 
iGrahamI ( 1200 lb . This at least confirms that our decomposi- 
tion method is sound, and able to recover the well established 
correlation. 

Pseudobulges and classical bulges overlap in range of bulge 
prominence, shown in the left panel of Fig. |6] The average 
classical bulge has {B/T) = 0.41. The width of the distribu- 
tion is 0.22. The average pseudobulge has {B/T) = 0.16. The 
width of the distribution is 0.11. Pseudobulges are on aver- 
age smaller fractions of the total luminosity of the galaxy in 
which they reside, ye t there is significant overlap (see also 
iDrorv & Fished |2007|) . This is no surprise, we would expect 
pseudobulges to be smaller as they are thought to be products 
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of secular evolution of the disk; hence their luminosity is ex- 
pected to be limited to some fraction of the luminosity of the 
disk. On the other hand, classical bulges form independently 
of their disk, presumably through major mergers. The full 
range values of B/T is therefore possible. 

NGC 2950 is a pseudobulge with B/T = 0.51. It is hard 
to believe that secular evolution can drive half of the stellar 
mass into the central region of a galaxy. This galaxy is an 
SBO wit h a prominent nu c lear b ar and nuclear ring, also no- 
ticed by lErwin & Sparkd (l2003l) . Thus, by our definition it 
is a good pseudobulge candidate. The next most prominent 
pseudobulge has B/T = 0.35; NGC 2950 is more than one 
standard deviation away in the distribution of pseudobulge 
B/T from the next most prominent pseudobulge. We feel that 
NGC 2950 is an exceptional case, and should not be taken 
as a normal pseudobulge. That it is an SO galaxy strengthens 
this interpretation. The unusually large B/T may be a result 
of the processes which made the galaxy SO (for e xample gas 
stripping by ram pressure and/or harassment; see Moo re etakl 
[1996) . rather than secular evolution. Also, our analysis may 
be an over simplification of the population of bulges, in that 
composite systems ("bulges" with both a pseudobulge and a 
classical bulge) may exists. This could artificially increase the 
B /T ratio. 

ICourteau et al.l (1 19961) show that the ratio of scale lengths 
for galaxies when fitted with a double expon e ntial ( bulge 
and disk) is tightly coupled. iMacArthur et al.l (l2003b find, 
with fits to Sersic plus outer exponential profiles in galax- 
ies spanning Hubble type Sab to Sd, that the coupling is 
(fe/h) = 0.22 ±0.09, and that the ratio becomes smaller to- 
ward later Hubble types (see also lower panel of Fig. |5]l. It 
appears that the correlation of r,,/!! with Hubble type may be 
driven primarily by the number of classical bulges in each 
Hubble type; this statement is very uncertain due to small 
numbers in each Hubble type, though. 

The right panel of Fig. |6] shows the distribution of r^/h 
for the galaxies in our sample (red lines are for classical 
bulges and blue lines represent the distribution of pseudob- 



ulges). We find that {r,/h) = 0.28 ± 0.28 for all galaxies 
in our sample, and that rg/h decreases toward later Hubble 
types. The avera ge is higher and the scatter is larger than in 
IMacArthur et al.' (2003), most likely because our sample tar- 
gets earlier-type galaxies. Considering bulge types separately, 
we find that for pseudobulges (rg/h) = 0.21 ±0.10, and for 
classical bulges (re/h) = 0.45 ±0.28. The distribution of rg/h 
for pseudobulges is clearly much narrower than that of the 
classical bulges. Furthermore, and h do not appear to be 
correlated in classical bulges. 

Our finding that the half-light radius of pseudobulges is 
well correlated with the scale length of the outer disk, while 
the scale length of classical bulges is not correlated with that 
of the disk is consistent with the interpretation that r^, cx /i is 
due to a secular formation of pseudobulges. However, inspec- 
tion of Fig. |6] shows that the range of rg/h over which we 
find classical bulges is large and overlaps with the pseudob- 
ulges significantly. Thus, rg/h may not be a good diagnostic 
tool for identifying pseudobulges and classical bulges; how- 
ever the finding that the scale length of the disk is correlated 
with the size of its pseudobulge but not correlated with clas- 
sical bulges is a strong indication that there indeed is a physi- 
cal difference between the formation mechanisms of different 
bulge types. 

3.2. The Sersic Index In Bulges 

Figure |2] shows the correlation of bulge properties with the 
shape of the surface brightness profile measured by the Sersic 
index, «i, of those bulges. We find that no classical bulge 
has Sersic index less than two, and very few pseudobulges 
(~ 10%) have Sersic index greater than 2. The correlations of 
Sersic index with other structural properties (r^, /ig, and M,,) 
appear to be different for pseudobulges and classical bulges. 
This suggests that the Sersic index may be used to distinguish 
pseudobulges from classical bulges. 

The top panel of Fig. [7] shows the correlation of Sersic in- 
dex with absolute V-band bulge magnitude for pseudobulges, 
classical bulges, and elliptical galaxies. Classical bulges fit 
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Fig. 7. — Correlation of the absolute magnitude (Mv), half-light radius 
(re), and suii'ace brightness at the half-light radius (pe) of bulges with bulge 
Sersic index. Pseudobulges are represented as blue crosses, classical bulges 
as red circles. We also show elliptical galaxies from Kormendy et al. (2006) 
for comparison, represented by black filled circles. The average uncertainty 
of the parameters of all bulges is represented by the en'or bars in the bottom 
right comer of each panel. 



well within the correlation set by elliptical galaxies. They 
cover a range of Sersic index from ^ 2-4.5, and their ab- 
solute magnitude covers the range My ^ - 1 7 .5 to -2 1 . Notice 
that the Sersic index of classical bulges is normally smaller 
than n = 4, suggesting t heir similarity to t he low -luminosity 
ellipticals discussed by iKormendv et al.l (l2006l) . Classical 
bulges, yet again, resemble little elliptical galaxies that hap- 
pen to be surrounded by a disk. The mean pseudobulge abso- 
lute magnitude is ^ 1 .2 magnitudes fainter than the average 
classical bulge magnitude, and the faintest bulges are all pseu- 
dobulges. Yet, pseudobulges do no have to be faint: classical 
bulges are not much brighter than the brightest pseudobulges. 
Our data cannot show if the brightest pseudobulges are bright 
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Fig. 8. — Residuals of pseudobulges (crosses), classical bulges (open cir- 
cles), and elliptical galaxies (closed circles) from a fit to the elliptical galaxies 
given by logC^f/kpc) = 2.91og(«ft) — 1.2. 

because of the nature of their stellar population, or because 
of greater stellar mass. Nonetheless, pseudobulges are not 
merely a change in observed properties at low magnitude. Fi- 
nally, we note that the overall shape of the My - n correlation 
for the super-set of systems represented in Fig.|7]is sim ilar to 
those f ound in previous st udies, [Grah am 2001; Balce lls et alJ 
2003; de Jong et alJl2004 . 

The bottom two panels of Fig. [T] show the correlation of 
Sersic index with the half-light quantities (radius containing 
half the light, r^, and surface brightness at that radius, /if) of 
pseudobulges, classical, bulges, and elliptical galaxies. The 
results are similar to those in the My - n plane. Classical 
bulges exist within the bounds set by elliptical galaxies, more 
specifically by low-luminosity ellipticals. Pseudobulges, on 
the other hand, populate the small n extreme of this correla- 
tion. The smallest objects, in radius, are pseudobulges. Yet, 
pseudobulges in general are not that much smaller on aver- 
age than classical bulges and low-luminosity ellipticals: the 
mean effective radius of pseudobulges lies within the range 
of effective radii spanned by the classical bulges and elliptical 
galaxies. 

To more strongly emphasize the difference between pseu- 
dobulges from classical bulges and ellipticals (hot stellar sys- 
tems), we show the residuals of all galaxies to the lCaon et alj 
(1 19941) relation r^, oc log(n/,) in Fig. [8] We fit this correla- 
tion to the elliptical galaxies only, and the resulting relation is 
log(r(,/kpc) = 2.91og(«i,)- 1.2. E-type galaxies and classical 
bulges have similar amounts of scatter Pseudobulges show a 
marked difference and systematically deviate further from the 
correlation as rih gets smaller whereas the scatter among clas- 
sical bulges is contained within the region occupied by ellip- 
ticals. Thus it appears that pseudobulges do not relate «/, to r^. 
in the same way as elliptical galaxies and classical bulges do. 
It is not clear from the bottom panel of Fig. |7] whether Sersic 
index correlates at all with effective radius for pseudobulges. 

Inspection of Fig. [T] shows that pseudobulges and classi- 
cal bulges have separate distributions of bulge Sersic index; 
the value that distinguishes the two types of bulges appears 
to be Hh ^ 2. Pseudobulges have Sersic index smaller than 
rth 2; classical bulges have larger values than rtt, ^ 2. We 
show this explicitly in Fig. |9] which shows the distributions 
of log(«i) for pseudobulges (blue lines) and the superset of 
classical bulges and elliptical galaxies (red lines), both binned 
to Slogirih) = 0.03. We also show a histogram of all ob- 
jects (pseudobulges, classical bulges, and elliptical galaxies) 
counted together. Each histogram is normalized by the num- 
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Fig. 9. — Histograms of bulge Sersic index, log(n;,), by bulge type. The 
top panel shows the histogram of all objects. The bottom panel shows pseu- 
dobulges (blue lines with mode log(M;,) ~ 0.22) and classical bulges + ellip- 
tical galaxies (red lines, mode log(ni,) ~ 0.54). We also overplot Gaussian 
distributions that fit the histograms. Histograms are normalized by the total 
number of galaxies in each group, so that each bin reflects the frequency of 
Sersic index within that object class. 



ber of objects in that group, so that the counts reflect the fre- 
quency of Sersic index within each class. Before discussing 
these distributions any further, we must strongly qualify this 
result. The sample studied in this paper is in no way complete. 
Thus, the relative abundance of pseudobulges and classical 
bulges cannot be determined from this distribution. However, 
we wish to study the distribution of Sersic indices in differ- 
ent bulge types. We do the best we can with the data we 
have at present; when relevant, we attempt to point out how 
our sample selection may bias results. We note that distribu- 
tions of bulge Sersic indic es in volume limi ted samples have 
been shown to be bimodal dAUen et alj20 06). We present any 
numbers derived from analysis of these distributions as pre- 
liminary. 

The top panel of Fig.|9]shows that the distribution of Sersic 
indices are clearly bimodal. This bimodality is completely 
coincident with the dichotomy in bulge morphology. The 
average values of bulges Sersic index for pseudobulges and 
classical bulges are 1 .69 ± 0.59 and 3.49 ± 0.60, respectively, 
and the average of the superset of classical bulges and el- 
liptical galaxies is 3.78 ± 1.5. The range of Hubble types 
chosen for this sample (SO-Sc) certainly biases the distribu- 
tion of Sersic index to higher values since the Sers ic index 
has been found to get smalle r at later Hubble types dGrahaml 
120011: iMacArthur et al.]|2003h . It is worth noting that the av- 
erage uncertainty of a bulge (pseudo- or classical) is compa- 
rable to the standard deviation of either of these distributions, 
{Anb) = 0.60. Inspection of Fig. |2]shows that neither log(rp) 
nor fie show a bimodal distribution. 

We fit a Gaussian to each histogram in Fig. |9] and solve for 
the Sersic index where the frequencies are equal. Because our 
sample is not volume limited, we weight the distributions such 
that the number of galaxies earlier than Hubble type Sc with 
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Fig. 10. — Two photometric projections of the fundamental plane. All 
symbols are the same as figure|2] 



pseudobulges is equal to the number of galaxies with classical 
bulges. We also solve for the case with 1/3 pseudobulges, and 
finally 2/3 pseudobulges. We find the critical Sersic index to 
be ticrii = 2.2 ±0.1. 

Of the 52 pseudobulges, 5 have Sersic index significantly 
above two (two pseudobulges have ^ 2). The interested 
reader can inspect their surface brightness profile in the ap- 
pendix, which shows fits to all surface brightness profiles. We 
take a few lines to discuss these galaxies here. NGC 4314 is 
an SBa galaxy with ni, = 2.37 ±0.78; this fit appears good, al- 
though a significant amount of the profile is not included due 
to the outer bar and nuclear ring. NGC 4258 is an oval galaxy, 
with nt, = 2.69 ±0.48. The outer oval ring affects the surface 
brightness profile of the outer disk, and thus the fit covers a 
narrow range. However, this is unlikely to affect the fit of 
the bulge too much, especially given the relatively low uncer- 
tainty in Sersic index in this fit. NGC 3627 is an Sb galaxy 
with tih = 2.90 ±0.83. Its bar is not easily detectable, and thus 
not removed from the fit, as we try not to remove any unneces- 
sary points. NGC 3642 is an Sb galaxy with = 3.37 ±0.61, 
the profile does not go very deep, and may be allowing nt to 
be artificially high. Finally, NGC 4826 has the largest Sersic 
index of any bulge classified as a pseudobulge, and the fit ap- 
pears quite good, unlike the others. We take these galaxies as 
exceptions rather than the rule, future studies comparing al- 
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ternative methods of pseudobulge detection (e.g. kinematics) 
will help shed light on their true nature. We remind the reader 
again that we are not classifying systems as composites, and 
these may have an effect. If for example a pseudobulge is em- 
bedded in a classical bulge (or vice versa), we might visually 
identify the pseudobulge even though the stellar light distri- 
bution is set by the classical bulge. Still, the dividing line 
of ni, ~ 2 is 90% successful at identifying pseudobulges, and 
thus is a good detector of pseudobulges. 

3.3. Fundamental Plane Projections 

Figure [To] shows the correlations of /ig - r^, and My - /ie, two 
proje c tions of the fundamental plane (iDiorgovski & DavisI 
119871; lDressleretal.lll987tlFaber etal.l ll989l). The top cor- 
relation shows the magnitude versus surface bright ness plane, 
and t he bottom panel shows the /i^ - relation (Kor mendyl 
[19771) . While the My versus plot has a lot of scatter, espe- 
cially in the pseudobulges, the radius-surface brightness plane 
shows significantly less scatter Yet in both of these funda- 
mental plane projections, those bulges which are further from 
the correlations established by the elliptical galaxies are pseu- 
dobulges. For those pseudobulges that deviate from the fun- 
damental plane correlation in Fig.[TO]that deviation is towards 
lower densities. 

This behavior has been noticed in the past. ICaroUol (Il999h 
finds that exponential bulges are systematically lower in effec- 
tive surface brightness than those better fit by an r'/^^-profile. 
We note that this could be an effect of only fitting either r'/"* 
or exponentials rather than a Sersic function to the bulge pro- 
files. If a bulge is not completely exponential then it may 
force other parameters to compensate f or the more restrictive 

S iararn eterization of the profile shape. iFalcon-Barroso et"a]l 
20021) find that bulges which deviate from the edge-on pro- 
jection of the fundamental plane are found in late type (Sbc) 
galaxies. However, at least with the galaxies in our sample, 
we cannot say unambiguously that this is a function of dif- 
fering bulge formation, because those bulges in Fig. [TOlthat 
deviate significantly from the correlations defined by the el- 
liptical galaxies are also very low B/T . Thus, it could be that 
the potential of the bulge is more affected by the outer disk. 



In fact, when looking at the pseudobulges in Fig. [TO] alone, 
one would not infer the presence of a strong correlation of 
either magnitude or effective radius with surface brightness. 
This is another indication that pseudobulges and classical 
bulges are different classes of objects. 

4. FLATTENING OF CLASSICAL BULGES AND 
PSEUDOBULGES 

There has been very little work on the distribution of flat- 
tenings of bulges albeit the structures present in pseudobulges 
(e.g. nuclear spirals or nuclear bars) suggest that pseudob- 
ulges should have higher angular momentum and thus be flat 
stellar systems. Howe ver, the end prod ucts of secular evo- 
lution need not be flat dKormendv'' 19931) . Bar buckling and 
unstable disks can both driv e stars higher above the plane of 
the disks (Pfenn iger & Norman 1990; Friedli 1999), thus cre- 
ating a pseudobulge that is les s flat than its associated outer 
disk. Nonetheless, the data of iKentI ( Il985h show that many 
bulges have median flattenings that are similar or greater than 
the median of the outer disk, and that flat bul ges ar e more 
common in late-type galaxies. iFathi & Peletie j te003b find a 
similar result that etuige/ ^disk > 0.9 in 36% of SO-Sb galaxies 
and 51% Sbc-Sm galaxie s (whe re e = 1 -b/a). However, note 
that iMQllenhoff & Heidi dioOl *) do not find this result, they 
find very few bulges are as flat as disks in ground based JHK 
imaging; 

Kormendv & KennicuttI ( 120041) include the flattening of 
bulges (as manifest through the ratio of bulge ellipticity to that 
of the disk) in their list of preliminary criteria for identifying 
pseudobulges. We can test this hypothesis with our sample. 
Figure [TT] shows the ratio of mean ellipticity of the bulge to 
that of the outer disk {{et,) / (e^)) for galaxies in this study. We 
do not include galaxies with average disk ellipticity less than 
0.2, as face-on projections of galaxies do not allow the flat- 
tening to be determined. Also, galaxies with B/T > 0.5 are 
removed; if the bulge dominates the entire potential it may 
set the shape of the disk and therefore affect the disk's ellip- 
ticity making it more like its own. To calculate the average 
ellipticity of the profile we use only the data points that are 
also included in the fit (filled circles in Figs. |4]and in the Ap- 
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pendix). It is a matter of interpretation as to what features are 
a part of pseudobulges, especially in light of the fact that some 
nuclei may be formed secularly along with the pseudobulge. 
However, as a matter of consistency we choose to focus only 
on those isophotes we call the "bulge" from the bulge-disk 
decompositions. The boundary between the bulge and disk 
is chosen as the radius at which the surface brightness of the 
Sersic function equals the surface brightness of the exponen- 
tial disk in the decomposition. Yet, in many galaxies contam- 
ination from the disk artificially raises the average ellipticity 
of the bulge. This contamination is evident in the ellipticity 
profile. Thus, we choose to average the bulge over a region 
in which there is little contamination present. Those radii are 
given in Table 2. 

The leftmost panel in Fig. [TT]compares (e/,) / (e,/) to bulge 
Sersic index for pseudobulges and classical bulges. It is quite 
evident that classical bulges, in our sample, are not as flat 
as pseudobulges. The average ratio of ellipticities of pseu- 
dobulges is 0.79 ±0.1, whereas the average ratio for classi- 
cal bulges is 0.49 ±0.14. In fact, the flattest classical bulge 
((e/,) / (ed) = 0.75) is less flat than the average pseudobulge. 

The middle and rights panels of Fig. [TT]compare this ratio 
to the bulge-to-total ratio and Hubble type, respectively. Clas- 
sical bulges show no obvious trend of (e/,) / (e^) with B/T or 



Hubble type. They are tightly clustered about the mean value 
~ 0.6 ±0.1. Pseudobulges however show a slight correlation 
of (eh) / (ed) with B/T. This correlation indicates that more 
prominent (in B/T) pseudobulges are slightly flatter. How- 
ever, the trend is we ak, and has large scat ter. Finally, we 
find a similar result to iKormendvl (1 1 9931) and lFathi & PeletieJ 
(|2003|) . namely that flatter bulges are more frequent in later- 
type galaxies. Furthermore, this trend only exists for pseu- 
dobulges. 

Invariably the flattest bulges are pseudobulges. However, 
we do not wish to overstate this result. It is possible that 
asymmetric dust extinction leads to higher apparent flatten- 
ing in the V band. Future work at near-to-mid IR wavelengths 
is likely to provide less ambiguous results. 

5. PSEUDOBULGES AS A CLASS OF OBJECTS 

In this paper we identify pseudobulges as having any, but 
not necessarily all, of several structures (nuclear bars, rings, 
and ovals). We treat bulges with these structures as a group. 
The motivation for doing so is that all these phenomena are 
similar to properties that are commonly associated with high 
specific angular momentum systems. At least in the sense 
that they are the complement of classical bulges they can be 
treated as a group. However, here we have to ask the follow- 
ing question: do thusly identified pseudobulges act as a single 
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Fig. 13. — Stnictural parameters of bulges as a function of the environ- 
ment. Field glaxies with classical bulges are marked by open squares, field 
pseudobulges with open diamonds. Cluster galaxies with classical bulges are 
marked by X's, cluster pseudobulge galaxies by crosses. All classical bulges 
are marked by red symbols, pseudobulges by blue symbols. 



class of objects or do significant differences exist among the 
objects we are identifying as pseudobulges? 

Figure[T2]suggests an answer to this question. The left pan- 
els in the figure show the distribution of bulge properties for 
bulges that have smooth isophotes (C=classical), bulges with 
nuclear spirals (nS), bulges with nuclear bars (nB), bulges 
with nuclear rings (nR), and bulges with a chaotic nuclear 
patchiness that resembles late type galaxies (nP). See Fig. |2] 
for example images of each of these features. In Sersic index, 
ellipticity ratio (bulge to disk), and bulge-to-total ratio there 
is no significant difference among the types of morphologies 
we call a pseudobulge. 

Bulges with chaotic nuclear-disk-like patchiness (nP) seem 
to have smaller Sersic index, and are more flat. It is unclear 
why this might be. This may be driven by them having smaller 
B/T and thus them being more affected by the disk potential. 
We also note that pseudobulges with «/, > 2 are almost all in 
bulges with nuclear spirals. It is possible that these bulges 
are not truly a spiral, instead some other phenomenon, like 
contamination of the light by the outer disk, is causing us to 
identify them as pseudobulges. However, in each parameter 
the distinction appears to be between classical bulges and the 
rest of the objects, rather than among the objects we call pseu- 
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dobulges. 

We also show Sersic index, ellipticity ratio (bulge to disk), 
and bulge-to-total ratio as a function of secular driving mech- 
anism. We separate galaxies as having no driving mecha- 
nism, an oval disk, or a bar Notice once again that we do 
not set grand design spiral as a class of object. There seems 
to be little differences in the averages, except ovaled galaxies 
have smaller Sersic index and flatter bulges. The phenomenon 
of flat bulges does not appear to require the presence of a 
bar. Thus it is not likely that all flat bulges are small bulges 
stretched by a bar potential making them flat. Also notice that 
the distributions of parameters for barred, oval, and normal 
spirals is roughly the same. Thus it does not seem likely that 
our method of removing the bar is artificially changing the 
bulge-disk decomposition parameters. 

6. ENVIRONMENTAL EFFECTS 

Secular evolution is not the only theory for building bulges 
that look like disks. Other possibilities such as extremely 
gas rich accretion events, distant gravitational encounters, or 
gravitational interactions with a cluster potential could all 
dri ve gas to the center of gal ax ies to increase B/ T, as f ound 
by iMastropietro et alj (l2005h . iKannappan et all (|2004 find 
that blue bulges are statistically more likely to have compan- 
ion galaxies within 100 kpc. Further, counter-rotating gas is 
frequently observed in spiral galaxies, and is taken as a sign of 
the galaxy having accreted galaxies in their past. Yet, pseu- 
dobulges are much more common in late-typ e galaxies, and 
the well known morphology-density relation dPressleiil 19801) 
shows that late-type galaxies are not found in dense environ- 
ments. Is it possible that pseudobulges in the field are formed 
through different mechanisms than those in cluster environ- 
ments? We can not know the answer to this for certain, how- 
ever we can look for signatures for such effects in our data. 

In Fig. [13] we show our primary result again, namely that 
pseudobulges have a lower Sersic index than classical bulges, 
but here we also mark environment. Field glaxies with clas- 
sical bulges are marked by open squares, field pseudobulges 
with open diamonds. Cluster galaxies with classical bulges 
are marked by X's, cluster pseudobulge galaxies by crosses. 
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There seems to be no significant differences in the structure of 
the surface brightness profile between the pseudobulges that 
reside in the field galaxies and the ones in cluster galaxies. 
Cluster pseudobulges are not higher mass nor systematically 
different in radial size than field pseudobulges. Nor are their 
surface brightness profiles preferentially steeper in clusters. 

The same holds for B/T: there is no substantial difference 
between pseudobulges in the field versus clustered environ- 
ments. If pseudobulge formation was driven primarily by 
tidal encounters with distant galaxies, one would expect that 
this effect should be enhanced in cluster environments, where 
such encounters are more frequent. This would result in more 
massive pseudobulges existing in cluster environments. Our 
sample does not seem to indicate that pseudobulges are any 
more luminous in cluster environments than in the field. We 
feel that this supports the notion that externally driven disk 
evolution is not likely the dominant affect in driving pseudob- 
ulge formation. 

Our results in no way rule out the possibility that pseudob- 
ulges are formed by gas-rich minor mergers. If we take an 
example of merging our Galaxy with one of the Magellanic 
clouds, and if this is done '--^1 Gyr, when the gas fractions were 
much higher, it is entirely plausible that the result could look 
similar to what we call a pseudobulge. However, it is not cer- 
tain if such a system would still be actively forming stars to- 
day, or if pseudobulge formation through minor mergers could 
only happen in extreme cases (e.g. prograde collisions at low 
inclination). Also, it is not clear how such accretion-driven 
formation of pseudobulges could maintain the bulge-disk cor- 
relations discussed in the previous sections. 

7. SUMMARY & DISCUSSION 

The main result of this paper is that bulge Sersic index, iib, 
is bimodally distributed in intermediate type galaxies where 
both classical bulges and pseudobulges exist. A value of 
~ 2 marks the boundary for separating morphologically- 
identified pseudobulges from classical bulges. Below tit, = 2 
no classical bulges are found, and above it very few pseudob- 
ulges are found. We also find that galaxies which are identi- 
fied as pseudobulges, using either bulge morphology or bulge 
Sersic index, are flatter than classical bulges. Thus suggesting 
that on average these systems are more disk-like in both their 
morphology and shape than are classical bulges. Pseudob- 
ulges exist that are as round as classical bulges, yet, invariably 
the flattest bulges are all pseudobulges. In both Sersic index 
and flattening ratio, our results show a homogeneity in classi- 
cal bulge properties, and a greater dispersion in pseudobulge 
properties. That is to say, classical bulges, in our sample, do 
not have Sersic index less than two, nor do they have high 
ratios of bulge-to-disk ellipticity. Pseudobulges, on average, 
have lower B/r than classical bulges, however we find pseu- 
dobulges with B/T extending to ^ 0.35. 

We find that the half-light radius of pseudobulges is well 
correlated with the scale length of the outer disk, while the 
scale length of classical bulges is not correlated with that of 
the disk. This is consistent with the interpretation that oc h 
is due to a secular formation of pseudobulges. Also, the fact 
that the scale length of the disk is correlated with the size of 
its pseudobulge but not correlated with the size of classical 
bulges is a strong indication that there indeed is a physical 
difference between the formation mechanisms of these differ- 
ent bulge types. 

In photometric fundamental plane projections, pseudob- 
ulges populate and extend the low luminosity end of the range 



occupied by classical bulges and elliptical galaxies. Pseu- 
dobulges deviate more from the relations set by the elliptical 
galaxies than classical bulges, and preferably towards lower 
density. In fact, it does not seem that on their own pseu- 
dobulges would establish any of these correlations, especially 

Mv-^Jie. 

In all correlations investigated in this work we find that the 
fundamental distinction is between classical bulges and pseu- 
dobulges. We do not find significant differences within the 
class of morphologies we identify as pseudobulges (nuclear 
spirals, nuclear bars, nuclear rings, and chaotic nuclear patch- 
iness). 

Is it possible that identifying pseudobulges visually using 
nuclear morphology is subject to an inherit flaw? That being 
that these systems are by definition bulge-disk systems, and 
thus is it possible that classical bulges coexist over large ra- 
dius with the central parts of a disk. Pseudobulges in SO galax- 
ies exemplify this concern. In many ways elliptical galax- 
ies and SOs are thought to be the extremes of a continuum 
of properties. Embe dded disks are known to exist in ellipti- 
cal galaxies (see, e.g. lScorza & Bendeiill995h . Would we call 
such a system where the embedded disk had a nuclear bar a 
pseudobulge? Most likely not. However, in SOs, which by 
definition have lower B/T, we are inclined to do so. Clearly, 
if there is no sharp change in properties along the E-SO contin- 
uum, there should be some composite systems with secularly 
formed structure in the disk and a classical (i.e. kinematically 
hot) bulge. (Note, though, that many SOs resemble defunct 
later type disk galaxies much mo re than they do resem ble el- 
liptical galaxies as suggested by Ivan den Berghlll976h . This 
may be the case in NGC 2950, yet it does not appear to be the 
rule. lErwin et"aLl ( |2004|) show that many bulges have more 
complicated dynamical profiles, with both hot an cold com- 
ponents. Yet, pseudobulges show a remarkable similarity to 
disk stellar populations, ISM, and star formation rates. These 
similarities ha ve been shown in the flatn ess of stellar popula- 
tion gradients (iPeletier & Balcelislll996h. and th e similarity of 
CO profiles to optical light (Regan et al.ll2001l) . iHelfer et all 
(2003) show that many CO gas maps of bulge-disk galax- 
ies have holes in their center, but also many CO profiles rise 
steeply to their center. A quick comparison of with our sam- 
ple shows th at those holes are found in classical bulges (e.g. 
NGC 2841). iFished (l2006f ) shows that the nuclear morphol- 
ogy of galaxies (as used here) predicts the shape of the 3.6-8.0 
micron color profile: disk-like bulges have flat color profiles, 
and E-like bulges have holes in 8.0 micron emission. If we 
were merely mistakenly identifying disks superimposed on 
bulges as physically different pseudobulges, we'd expect to 
find larger B/T for such systems as compared to the bulges 
where we do not see evidence for the presence of a disk which 
we call classical in this work. This is because of the added 
light of the disk onto a fraction of bulges drawn from the same 
underlying distribution. However, we do not find this to be the 
case. Fig. |6] shows that pseudobulges are on average smaller 
than classical bulges. 

It is not clear whether pseudobulges participate in the var- 
ious correlations presented in this paper. If we look at the 
pseudobulges alone, they only show a convincing correlation 
in the Kormendy relation (/i^, - ). Their distribution better 
resembles scatter diagrams in the My-fi^, /^e-n, and r^-n 
planes. However, it would seem to be an odd coincidence that 
in every parameter combination pseudobulges just happen to 
fall in the right location to extend the correlations set by the 
classical bulges and the elliptical galaxies. Putting all this 
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Fig . 15. — Bulge-to-total ratio (lower panel) and ratio of the bulge half-light 
radius to the scale length of the disk (top panel) plotted against the logarithm 
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together implies that we do not really understand the details 
of pseudobulge formation very well. It is also possible that 
the decomposition of pseudobulges may not be an appropriate 
procedure. Bulge-disk decomposition assumes that the struc- 
tures are distinct (and that light extends to radius of infinity). 
If we accept that pseudobulges are formed through disk phe- 
nomena, then such a bulge-disk decomposition may not be an 
adequate description of those systems. In reality they might 
simply be a component of galactic disk. A lso, parameters de- 
rived from the Sersic function are coupled (iGraham & CoUessI 
[l997). Thus, treating pseudobulges in the same way as clas- 
sical bulges may artificially force them to extend some photo- 
metric correlations. 

Concerning the six bulges with disk-like morphology near 
their centers that have rib > 2, two have (et,) / {ej) > 0.75 
(the maximum of classical bulges), and NGC 4314 has a 
very prominent star forming nuclear ring, making this bulge a 
strong candidate for being a pseudobulge. This galaxy might 
be an example for a composite system. Their existence cer- 
tainly underscores the value of having as much information 
as possible when diagnosing bulge types. Future work may 
prove illuminating. For example, to what extent should we 
trust bulge morphology as an indicator of secular evolution? 
It is possible for the human eye to mistake merely the pres- 
ence of dust for spiral structure? 

We also find that the ratio of scale lengths, re/h, in pseu- 



dobulges is more tightly correlated and closer to thos e values 
reported by other authors (e.g. Mac Arthu r et alj|2003h than in 
classical bulges. However, the range in values of r^/h for clas- 
sical bulges and pseudobulges are similar; therefore this does 
not provide a good diagnostic tool for finding pseudobulges. 
Pseudobulges extend the parameter correlations of photomet- 
ric quantities (r^,, fig, My, and tih). However in many of these 
parameters it is unclear if pseudobulges actually show a cor- 
relation on their own. 

We can compare the structural properties of the bulges in 
our sample the output of simulation. Unfortunately simula- 
tions of galaxies that resemble real galaxies including stars, 
dark matter, gas and star formation (and possibly feedback) 
are quite difficult. Thus there does not exist a statistically rel- 
evant set of simulations for full comparison. Nonetheless, we 
c an compare our r e sults t o those that currently exist. 

iDebattista et al.l (|2004|) provides a set of simulations that 
generate bulges from bar buckling in pure stellar (no gas) sys- 
tems. The resulting bulges typically have nt, ^ 1.5, which 
is consistent with what we find. The bulges in our sample 
tend to be more round than their associated outer disk, and 
IDebattista et al.l (|2004|) separate their bulges based on flat- 
tening ratios. So we will only compare to those bulges that 
are less round than the outer disk. Their simulated galaxies 
have B/D = 0.2-1.0 (where B/D is bulge-to-disk Hght ra- 
tio). Where as our sample has a median B/D = 0.12 with a 
standard deviation of 0.20. Further, they are able to recover 
the coupling of r^, and h that we find in pseudobulges. Thus 
our pseudobulges tend to be a bit smaller. It may be that bar- 
buckling is one way to make a pseudobulge, as indicated by 
the the fact that those simulated bulges from IDebattista et all 
( |2004|) are contained within the set of pseudobulges, but do 
not span the whole range in properties. 

We reiterate a staternent by Andreda kis] (| 19981) . that is also 
discussed in iKormendv & Kennicutti (120041) . namelv that we 
do not really understand why pseudobulges should have a cer- 
tain Sersic index. It is understandable t hat mergers would 
driv e up the Sersi c index, as discussed by van Albada' dHH 
and iKormendv et al. (2006), and thus classical bulges are 
found with higher Sersic index. lAguerri et alJ (120011) simu- 
lates minor mergers (accretion of satellites on the order of the 
mass of the bulge) in bulge-disk galaxies. They show that 
the Sersic index grows as the amount of mass accreted be- 
comes larger Eliche-Moral et al. (2006) takes the study of 
sat ellite accretio n to lo wer densities than those simulations 
of lAguerri et al.l (1200 ll) . Their simulations do not include 
as an d star formation, and thus as with the IDebattista et alJ 
20041) simulations they should be read with that caveat in 
mind. They find that low density satellite accretion does 
not necessarily drive Sersic index above the critical value we 
fmd(ncrit ^ 2). Yet, they also find that satellite accretion leads 
to a simultaneous increase in B/T. 

In Fig. [15] we therefore plot the bulge-to-total ratio, B/T, 
and the ratio of the bulge half-light radius to the scale length 
of the disk, r^/h, against «/,. There is a tight correlation for 
classical bulges. However, pseudobulges do not follow this 
correlation. In fact, neither B/T nor rg/h correlate with «/, in 
pseudobulges at all in the range of parameters shown, span- 
ning an order of magnitude in B/T. The absence of correla- 
tion between rii, and B/T in pseudobulges seems to be point- 
ing to a non-merger driven formation scenario for pseudob- 
ulges. Yet, these results are suggestive at best. There may 
be some underlying correlation that is destroyed by other fac- 
tors (e.g. gas fraction or collision parameters). More work is 
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needed from both simulations and observations. That we do 
not see correlatio n of rib and B/ T for low-Sersic index bulges 
(as found in Eliche-Moral et al. 2006) may be indicating that 
the classical bulges we observe today are the products of mul- 
tiple mergers, and possibly at higher redshifts there is a pop- 
ulation of low-Sersic index classical bulges. Conversely it is 
entirely possible that a population of small classical bulges 
exists, yet they are embedded within pseudobulges. In this 
case the mass of a bulge, and hence the B/T within a specific 
galaxy would be coming from multiple mechanisms. 

If pseudobulges are built by secular evolution, the simplest 
mechanism controlling the Sersic index in bulges is that nb 
grows with time as the bulge-to-total ratio increases. How- 
ever, as discussed above, it appears that there is not a strong 
connection between B/T and Sersic index. This would im- 
ply that Sersic index in pseudobulges is not a time-dependent 
quantity. We wish to emphasize, though, that the error in the 
measurement of is large, and might be masking an underly- 
ing weak correlation. Also, by focusing on intermediate-type 
galaxies, where both pseudobulges and classical bulges oc- 
cur, we miss very low B/T systems. Expanding a sa mple t o 
later types might reveal a weak correlation (Graham 200 lb . 
Since r^, correlates tightly with li for pseudobulges, it is no 
surprise that r^/h does not correlate with nb in pseudobulges. 
The linea r coupling of and h is well established in late type 
galaxies (ICourteau et al.lll996h . 

It may be that the Sersic index of pseudobulges is merely 
another manifestation of the dynamical state of the system. 
Stars in bulges with larger amounts of random motion often 
take on radial orbits, thus climbing higher out of the potential 
well. Thus there is more light at large radius, increasing n^. 
However, stars in orbits with higher amounts of angular mo- 
mentum would be less likely to take on radial orbits and thus 
bulge light would contribute less at large radius; driving nb 
down. The observation that the distribution of Sersic indices 
in bulges of galaxies from E to Sc is bimodal then strength- 
ens the claim that what we are calling pseudobulges are not 
merely the low-mass counterparts of the same phenomena that 
form classical bulges and elliptical galaxies. This description 
fits well with the theory that secular evolution forms pseudob- 
ulges, and mergers, whether by a single event or succession 
of events, form classical bulges. The higher angular momen- 
tum (and thus low Sersic index) systems have not had major 
mergers and thus have not experienced the violent processes 
that lower the ratio of rotational velocity to velocity disper- 
sion. What is left unknown is why pseudobulges are not ex- 
ponential, and also what keeps them from having larger Sersic 
indices. 

We conclude by returning to the results of Fig. |6] and the 
implications of the range of B/T for pseudobulges. It is 
now well known that the presence of many low -B/T sys- 
tems pre sents a problem for current A CDM galaxy formation 
theories (D'On ghia & Burkertll2004l) . In our sample we find 
no classical bulge galaxies with B/T < 0.1. Recent studies 
which compare merger histories in ACDM simulations to the 
observed frequency of bulgeless galaxies suggest that either 
there are too many mergers in the simulations, or that disks 
must be much more robust to the merg ing process th an pre- 
viously thought (Stewart et al. 2007; Ko da et al.|[2007l) . Pseu- 
dobulges span a range of B/T from 0.35 to zero. If pseudob- 
ulges form through internal-disk processes then a galaxy with 
a pseudobulge can be thought of as a pure disk galaxy. Thus 
current estimates of the number of low -B/T systems could 



only be thought of as lower-limits; the existence of pseudob- 
ulges would make the problem of forming bulgeless systems 
even more pressing. 
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